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' ^ [ Abstract. We establish an analogue of the Bochner theorem for first order operators of 

["T | ■ Dunkl type, that is we classify all such operators having polynomial solutions. Under natural 

conditions it is seen that the only families of orthogonal polynomials in this category are 
' limits of little and big (j- Jacobi polynomials as q = — 1. 
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■ 1 Introduction 

> 

| We have introduced recently new families of "classical" polynomials \15\ [16] through limits of 

little and big q- Jacobi polynomials when q goes to —1. The term "classical" is taken to mean 
that the polynomials P n (x) are eigenfunctions of some differential or difference operator L 

LP n (x) = X n P n (x). (1) 

In the case of the novel classes of polynomials that we have identified, the operator L is of first 
order in the derivative operator d x and contains moreover the reflection operator R defined by 
^ . Rf{x) = f(—x). In other words the operator L can then be said to be a first order differential- 

difference operator of Dunkl type [6]. 

The little —1- Jacobi polynomials depend on 2 real parameters a, (3 such that a > —1, B > —1. 
They are orthogonal on the interval [—1,1] with respect to the positive weight function [T5] 



w (x) = (x + l)(l -x 2 ) {a ^l 2 \xf. (2) 

The big — 1-Jacobi polynomials depend on 3 real parameters a, (3, c such that a > —1, (3 > — 1, 
< c < 1 and are orthogonal on the union of two distinct intervals [— 1, — c] and [c, 1] with 
respect to the positive weight function [16] 

w ( x ) = M( x + l)( x _ c ) (l _ x 2 )( a-l)/2/ 2 _ c 2)03-l)/2_ (3) 
X 

When c = the weight function ([3]) becomes ([2]). This corresponds to the well known limit 
relating the big to the little q- Jacobi polynomials [9]. 
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Both weight functions ([3]) and ([2]) are not associated to ordinary classical polynomials. Ne- 
vertheless, the weight function corresponding to the little — 1-Jacobi polynomials, belongs 
to the class of so-called generalized Jacobi weights (see, e.g. [12]). 

Like the ordinary classical orthogonal polynomials, the big and little —1- Jacobi polynomials 
possess many remarkable properties, in particular, they can be presented in terms of Gauss' 
hyper geometric functions [15| [T6] . 

The purpose of the present paper is to examine in general the eigenvalue problem ([T]) for 
first-order operators L of Dunkl type and hence to classify the families of polynomials that are 
eigenfunctions of such operators. As a result, under the natural condition that the operator L is 
potentially self-adjoint, we show that the big — 1-Jacobi polynomials (and their little — 1-Jacobi 
polynomials limit) exhaust the list of orthogonal polynomials with that property. 

Note that these polynomials have also been found to arise from the Bannai-Ito polynomials pQ 
when the orthogonality support in the latter case is extended to an infinite set of points in a limit 
process [16]. 

2 Dunkl-type operators of first order 
and their polynomial solutions 

Consider the most general form of linear differential operators L of first order which contain 
also the reflection operator R: 

L = F (x) + F x (x)R + G (x)8 x + Gi(x)d x R, (4) 

where Fq(x), F±(x), Gq(x), G\(x) are arbitrary functions. 

We are seeking polynomial eigensolutions of the operator L, i.e. for every n we assume that 
there exists a monic polynomial P n (x) = x n + 0(x n ~ 1 ) which is an eigenfunction of the opera- 
tor L with eigenvalue A n ([I]). In what follows we will assume that 

A n ^ for n = 1, 2, . . . , A n ^ A m for n ^ m. (5) 

We first establish the necessary conditions for the existence of such eigensolutions. 

Consider the monomials x n , n = 0,1,2,.... If condition ([TJ holds for all n, the opera- 
tor L must preserve the space of polynomials of any dimension, i.e. for any polynomial ir(x) 
of degree n, we have Ltt(x) = 7r(x), where 7r(x) is a polynomial of the same degree as ir(x). 
Indeed, since ir(x) can be expanded as a finite linear combination of the polynomials Pi(x), 
i = 0, 1, ... ,n, the condition X n ^ guarantees that the polynomial tt{x) = Lir(x) and ir(x) 
have the same degree. 

Hence from condition (pQ) it follows that 

Lx n = Q n (x), n = 0,1,2,..., (6) 
where Q n (x) are some polynomials of degree n: 

Q n {x) = \ n X n + 0{x n ~ 1 ). 

Conversely, assume that condition (|6|) is valid for all monomials x n , where Q = n{x) are some 
polynomials with leading term A n . 

Then it is easily seen that for any n = 0, 1,2, ... there exists a monic polynomial P n (x) = 
x n + 0{x n ) satisfying the eigenvalue equation ([1]). 

Indeed, from ([6]) and from conditions ([5]) for A n it follows that the operator L transforms 
any linear space of polynomials of exact degree n into itself. Hence, from elementary linear 
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algebra it is known that there always exists at least one eigenvector with the eigenvalue X n . 
This eigenvector corresponds to the polynomial solution P n (x) of given degree n. Moreover, 
again by conditions ([5]), it follows that for every n = 0, 1, 2, . . . such solution is unique. We thus 
proved that conditions ([T]) and © are equivalent (restrictions ([5]) are assumed, of course). 

We shall thus proceed to verify condition © for the first values of n, i.e. n = 0, 1, 2, 3. This 
will provide necessary conditions for the functions Fi(x), Gi{x) to satisfy. 

Consider the case n = 0. Without loss of generality we can assume that 

L{1} = 0. 

Indeed, by condition ([6]) we have £{1} = c for some constant c. Since this constant can be 
incorporated additively into the function Fq(x), we can assume that c = 0. Then, from © we 
find that necessarily F\(x) = —Fq(x), i.e. that the operator L must have the form 

L = F(x)(I -R) + G (x)d x + G 1 (x)d x R, (7) 

where F(x) = Fq(x) and / is the identity operator. In what follows we will use the form ((7|) for 
the operator L. 

Consider now the case n = 1. We have L{x} = Q\{x) with Qi(x) a linear function. Hence, 
we obtain the condition 

2xF(x) + G (x)-G 1 (x) = Q 1 (x). (8) 

Analogously, for n = 2, 3 we obtain the conditions 

2x(G (x) + G 1 (x)) = Q 2 (x) (9) 

and 

2x 3 F(x) + 3x 2 (G (x) - Gx(x)) = Q 3 (x). (10) 

We thus have 3 equations ([S]), Q, (fTU|) for the 3 unknown functions F(x), Go(x), G\{x). Solving 
these equations, we can present their general solution in the form 

„ , > u v o , . a v\ 

G (x) = ^H \-p + T x, &i(x) = 5-H \-pi+nx, 

^) = -4 + ^ + i +^ (id 

x A 2x z x 

where //, vq, v\, po, pi, tq, n, ^, r\ are arbitrary constants. In what follows we will assume the 
restrictions 

ti^±t , 2 V + (2N + l)(r - n) ^0, iV = 0,l,2,... (12) 
We thus have 

Theorem 1. Let p, v\, po, p\, tq, t\, ij be arbitrary complex parameters subjected to 
conditions (|12p . Construct the operator L using formulas ([7]) and (jlip . Then, for any n, 
condition ([6]) holds with 



X, 



(to + Ti)n if n is even, 

2?? + (to — T\)n if n is odd, 



and the operator L has a complete set of polynomial eigenf unctions LP n (x) = \ n P n (x), n 
0, 1, 2, . . . with the nondegeneracy conditions ([5]). 
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The proof of the theorem follows from the simple observation that for any n = 1, 2, 3, . . . we 
have 

Lx n = X n x n + 

(i) 

where the coefficients k\ , i = 1, 2, 3 are straightforwardly found from the explicit expression of 
the operator L. Hence the operator L preserves the space of polynomials. The nondegeneracy 
conditions © are equivalent to the conditions (fl~2|) . 



3 Symmetrizable operators L 

In the previous section we identified the most general form of first-order operators L of Dunkl- 
type that have nondegenerate polynomial eigensolutions P n (x). Assume now that the poly- 
nomials P n (x) are orthogonal, i.e. that there exists a nondegenerate linear functional a such 
that 

(a,P n (x)P m (x)) = h n 5 nm (13) 

with some nonzero normalization constants h n . 

This condition means that the operator L is symmetric with respect to the functional a: 

(a, L{V{x)}W{x)) = (a, L{W(x)}V(x)), (14) 

where V(x), W{x) are arbitrary polynomials. Property (|14[) trivially follows from ([1]) and the 
completeness of the system of polynomials P n (x). 

In the case of a positive definite functional, i.e. h n > 0, there is a realization of (|13p in terms 
of an integral 



{a,P n (x)P m (x)) = [ P n (x)P m (x)da(x) 

J a 



with respect to a measure da(x), where cr(x) is a nondecreasing function of bounded variation [8]. 
The limits a, b of integration may be finite or infinite. 

We restrict ourselves to the case where there is a positive weight function w(x) = da(x) /dx>0 
inside the interval [a, b] (we do not exclude the existence of concentrated masses at the end- 
points a, b of the interval). 

Under such restrictions we have the following necessary condition for the operator L: 

{w(x)L)* = w(x)L, (15) 

where M* stands for the Lagrange adjoint operator with respect to the operator M. Moreover, 
in what follows we will assume that all functions F(x), Gq(x) and Gi(x) are real for the real 
values of the argument x. Recall that if the operator M is a differential operator 

N 



k=0 

with real- valued functions Ak(x) then the Lagrange adjoint operator is defined by 

N 



M* = Y,(-l) k d k x A k (x) 

k=0 
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In the presence of the reflection operator R, i.e. in the case where the operator M has the form 

N N 

M = J2 A k (x)d k x + B k {x)d k x R 

fc=0 k=0 

with some real- valued functions A k (x), B k (x), we define the Lagrange adjoint operator by 

N N N N 

M* = Y,(-^ k d k x A k (x)+Y,(-^) k Rd k x B k (x) = J2(-^ k d k x Mx) + Yd k x B k (-x)R, 

k=0 k=0 k=0 k=0 

where we used the following formal rules 

Rd x = -d x R, R* = R, RB(x) = B{-x)R. 

These formulas assume that the interval of orthogonality is necessarily symmetric, i.e. b = —a 
(note that this includes possible situations where the interval of orthogonality is the union of 
several intervals; in such a case, pairs of corresponding intervals should be symmetric with 
respect to the reflection x — > —x). 

We thus restrict our considerations to operators L which possess property (|15p with some 
positive function w(x) inside the symmetric interval (or intervals) of orthogonality. Such ope- 
rators are called symmetrizable while the corresponding function w{x) is called the symmetry 
factor for L [IHE]. 

It is easy to see that this symmetrizability property excludes the existence of the term Go(x)d x 
in (|4|) (this statement is valid only for the case of the real valued functions F(x), Gq(x), Gi(x) 
which has been assumed). Hence we restrict ourselves in the following, to operators of the form 

L = F(x)(l -R) + Gi{x)d x R 

with Gq(x) = 0. 

This implies p = v$ = po = t$ = and from (llip . we obtain the following expressions for the 
functions G\{x) and F{x): 

Gi{x) = — + pi + nx, F(x) = -p^ + - + r] 
x 2x z x 

with 5 arbitrary parameters ui, pi, t±, £, r/. 

Note that we can multiply the operator L by an arbitrary nonzero constant L —¥ kqL. A sca- 
ling transformation of the independent variable x — > K\x with an arbitrary nonzero constant K\ 
can also be performed. Using this freedom, without loss of generality, we can always reduce the 
number of arbitrary parameters to 3. 

Assume that v\T\ ^ 0. We can then rewrite G\{x) in the form 

Gl (x) = gl(x - d)(x + c) , (16) 

x 

where gi, d, c are real constants. Under the condition that d ^ c, we can set g\ = 2 and d = 1 
by an appropriate choice of the parameters Ko, k\ so as to have 

GM = 2( *" 1)( * + C) (17) 

x 

with one arbitrary real parameter c. We will assume that < c < 1. (The cases c> 1 or c < - 1 
can be analyzed analogously.) For the function F(x) we have 

F(x) = -4 + ^-^- a -/3-l (18) 

x z x 

with two arbitrary real parameters a, (3, c. 
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In [16] it was shown that the operator L with the coefficients G\(x), F(x) given by (|16[) 
and (|18p have big — 1-Jacobi polynomials as eigensolutions. In the limiting case c = the little 
— 1-Jacobi polynomials are obtained. 

Hence we have 

Theorem 2. The functions G\{x) and F(x) given by (|16p and (|18p thus provide the most 
general first order operator of Dunkl-type under our hypotheses. They define the operator L that 
has the big —1-Jacobi polynomials as eigenf unctions 

Note that we have considered the generic (nondegenerate) choice of real parameters in the 
functions G\(x) and F{x). In the next section we will examine the symmetrizability property 
of the operator L and consider also the degenerate cases with respect to the parameters. 



4 Pearson-type equation for the weight function 

It is instructive to derive directly the expression for the weight function w(x) from condition (115p . 

Denote M = w(x)L = w(x)F(x)(l — R) + w{x)G\{x)d x R. Assuming that all parameters of 
the functions F(x), G\(x), w(x) are real, we have 

M* = w(x)F(x) - w(-x)F(-x)R - d x Rw(x)Gi(x) 

= w{x)F(x) - w(-x)F(-x)R + d ( w (- x \ Gl (- x » R + wt-x^GA-x^dxR. 

ax 

Condition M* = M thus means that the following two equations hold: 

w(x)G 1 (x) = w(-x)G 1 (-x) (19) 

and 

w(-x)F(-x) - w(x)F(x) = ^-w(-x)G 1 (-x) = —w(x)G 1 (x). (20) 

ax ax 

With G\(x) given by (JTTJ) , we have from ([!§]) 

w(x) = 9(x)(x + l)(x - c)W(x 2 ), (21) 

where 9{x) = x/\x\ is the standard sign function and W(x 2 ) is an arbitrary even function of x. 
We will suppose that W(x 2 ) > inside the interval [c, 1]. Assume that a > — 1, f3 > — 1. Then 
from (|2ip it follows that the function w{x) is positive inside the two symmetric intervals [—1, — c] 
and [c, 1] . 

Substituting expression ()21|) into the second condition of (|20p and using the explicit expres- 
sion (|18p . we obtain a differential equation for the function W(y): 

2(y -l)(y- c 2 )W'(y) + ((2 - a - p)y + /3 - 1 + c 2 (a - l))W(y) = 0, 
which has the following general solution 

W{y) = C{l-y)^- 1 ^ 2 {y-c 2 f~ l)l2 
with an arbitrary constant C. Thus the function w(x) has the expression 

w{x) = C9{x){x + l)(x - c)(l - x 2 )( a -V /2 (x 2 - c 2 ) (/3 ~ 1)/2 . 

Assuming C > and a > — 1, /3 > — 1, we see that w(x) > inside the two symmetric intervals 
[— 1, — c] and [c, 1]. This coincides with the weight function w(x) corresponding to the big 
— 1-Jacobi polynomials which was derived in [16] from a completely different approach. 
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Consider now possible degenerate cases of the function G\{x). 

(i) Assume that v\ = but p\T\ 7^ 0. Then using scaling transformations, we can always 
reduce the functions G\{x) and F{x) to 

Gi(x) = 2(1 -x), F(x) = a + p + 1 - p/x. 

From (|19p we get w(x ) = (x + l)W(x 2 ) with some function W(y). Then from (|2ip we obtain 
a Pearson-type differential equation or the function W(y): 

2y(y - l)W'(y) + ((1 - a - p)y + p)W(y) = 0. 
Its general solution is 

W(y) = Cyf"\l-y)W. 
Hence we obtain for the weight function 

w(x) =C{x + l)\x\> 3 (l-x 2 ) ia ~ 1)/2 . 

This weight function corresponds to the one of the little — 1-Jacobi polynomials (|2J) which are 
orthogonal on the interval [—1,1] |15] . Equivalently, this case corresponds to setting c = in 
formulas (JUJ) and (fl8]) . 

{ii) Assume that u± = pi = but ri 7^ 0. Then we can write without loss of generality 

G 1 {x) = 2x, F{x) = a + p + 1- p/x. 
From (I19j) we now have w{x) = 6{x)W(x 2 ) and from (120j) follows the Pearson equation 

2yW'(y) + (a + p + 2)W(y) = 
with solution = Cy-( Q+/3+2 )/ 2 . Whence 

to(z) = C9{x)\x\- {a+l3+2)/i . (22) 

The function w(x) of (|22p is not positive however, in any symmetric interval [—a, a] (or any 
union of symmetric intervals). Hence, this case does not lead to orthogonal polynomials with 
positive measure on the real axis. 

(Hi) Consider the degenerate case when the numerator in expression (I16p has coinciding zeros. 
Using scaling transformations, G\{x) and F{x) can be cast in the form 

2(r — 1 ) 2 

Gi(x) = -±- L , F{x) = x- 2 + ax~ 1 + b 

x 

with arbitrary constants a, b. 
From condition (1191) we find 



w(x) = 8(x)(x + l) 2 W(x 2 ) 

with an arbitrary even function W(x 2 ). From condition (|20j) we then obtain a differential 
equation for the function W{y): 

2(y - l) 2 W'{y) + ((b + % + 2a + b - 1) W(y) = 

with general solution 

W(y) = C(y-l)-^/ 2 e XV (^±± 

\ y- 1 
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with an arbitrary constant C. Thus the weight function can be presented in the form 

w(x) = C8(x)(x - I) 2 (x 2 - l)" (b+3)/2 exp fe^) ■ 

It is easily seen that due to presence of the sign function 9(x) there are no symmetric intervals 
[— d, d] (or pairs of symmetric intervals) inside which the function w(x) can be positive. Hence 
in this (degenerate) case the operator L is not symmetrizable. 

(iv) Consider the case T\ = 0, v\p\ ^ 0. We can here posit G\{x) = 2(1 — x _1 ) and 
F(x) = —x~ 2 — ax -1 — (3 — 1. This leads to w(x) = 6(x)(x + l)W(x 2 ) and for the function W(y) 
we have the equation 

2(1 -y)W'(y) + (a + (3)W(y) = 0, 

whence 

w{x)=C9(x)(x + l){l-x 2 ) {a+P)/2 . 

Again it is impossible to realize the condition w(x) > on a symmetric interval. 

(v) Finally consider the case when t\ = p\ = 0. ThenGi(x) = — 2/x, F(x) = — x~ 2 — ax -1 — /3. 
This case corresponds to the weight function 

w(x) = C9(x)exp ^— ^x 2 

and here also it is impossible to have a positive weight function on a symmetric interval. 

We see that in all cases, the function w(x) possesses a factorization w(x) = ir(x)W(x 2 ), where 
ir(x) is a polynomial of first or second degree, and the function W(y) obeys a Pearson- type equa- 
tion that we have derived. Roughly speaking, the function W(y) resembles the corresponding 
weight function for the classical orthogonal polynomials |13j . But, in contrast to the classical 
case, no analogs of Laguerre or Hermite polynomials appear in the case q = —1. 



5 Conclusions 

We showed that the symmetrizability property of the Dunkl operator L with a positive weight 
function w(x) on a symmetric interval (or a union of symmetric intervals) leads to very strong 
restrictions upon the associated "classical" orthogonal polynomials. Namely, only big and little 
— 1-Jacobi polynomials admit Dunkl analogs. A number of generalizations offer themselves. 

For example it was shown in [3] that the generalized Hermite polynomials H^' \x) and the 
generalized Gegenbauer polynomials Sn'^ (x) (see, e.g. [H [T^l [2] for their definition and pro- 
perties) are eigenfunctions 

L fl P n ((x) = \ n P n (x) 
of an operator which is quadratic 

L/j, = v(x)T 2 + r(x)T M (23) 
with respect to the Dunkl operator 

T^ = d x + fix~ 1 (I-R). 



Here a(x) is a polynomial of degree not exceeding 2 and deg(r(x)) = 1. 
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Hence the generalized Hermite and Gegenbauer polynomials can be considered as "classical" , 
but in contradistinction to the general situation considered here, they satisfy an eigenvalue 
equation of second order in the Dunkl operator. The classification of all second-order Dunkl- 
type operators having orthogonal polynomials as solutions is a more involved problem. 

Note that the operator L defined by (f23|) is a direct generalization of the classical hypergeo- 
metric operators 

L = a{x)d 2 x + r{x)d x 

leading to the classical orthogonal polynomials as eigensolutions [13j . 

Let us mention also that operators of Dunkl-type have also been used in [10] and [5] to 
construct some polynomial eigenvalue solutions. In these cases however, the Dunkl operators are 
not symmetrizable, and the solutions therefore do not belong to the class of ordinary orthogonal 
polynomials. 
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